Astronomy & Astrophysics manuscript no. ms 


©ESO 2010 


March 3, 2010 





Studying the small scale ISM structure with supernovae 

F. Patat^ N.L.J. Cox^-^, J. Parrent'^ ^ and D. Branch'^ 

' European Organization for Astronomical Research in the Southern Hemisphere (ESO), K. Schwarzschild-str 2, 85748, Garching 
b. Munchen, Germany 
e-mail: fpatat@eso.org 

- Herschel Science Centre, European Space Astronomy Centre, ESA, PO Box 78, 28691 Madrid, Spain 
^ Institute of Astronomy, K.U. Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium 

Department of Physics and Astronomy. University of Oklahoma, Norman, OK 73019 
^ Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755-3528 

Received 08/01/2010; accepted 25/02/2010 

ABSTRACT 

Aims. In this work we explore the possibility of using the fast expansion of a Type la supernova photosphere to detect extra-galactic 
ISM column density variations on spatial scales of ~100 AU on time scales of a few months. 

Methods. We constructed a simple model which describes the expansion of the photodisk and the elfects of a patchy interstellar cloud 
on the observed equivalent width of Nai D lines. Using this model we derived the behavior of the equivalent width as a function of 
time, spatial scale and amplitude of the column density fluctuations. 

Results. The calculations show that isolated, small (<100 AU) clouds with Nai column densities exceeding a few 10" cm"- would be 
easily detected. In contrast, the effects of a more realistic, patchy ISM become measurable in a fraction of cases, and for peak-to-peak 
variations larger than ~10'^ cm"- on a scale of 1000 AU. 

Conclusions. The proposed technique provides a unique way to probe the extra-galactic small scale structure, which is out of reach 
for any of the methods used so far. The same tool can also be applied to study the sub-AU Galactic ISM structure. 
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1. Introduction 

For many years it was accepted that the minimum size for the 
column density fluctuations in the Galactic ISM is around 1 pc 
(~2xl0^ AU). The common understanding was that although 
sub-parsec structures do exist, only a tiny fraction of the col- 
umn density could be ascribed to these small scales (Dickey 
& Lockman J 9900 . However, the pioneering VLBI work by 
Dieter, Welch & Romney ( 119761) . and the later confirmation by 
Diamond et al. ( 119891 ) demonstrated the existence of significant 
fluctuations over scales of 20 AU. These findings were con- 
firmed by 21 cm absorption measurements against high-velocity 
pulsars (Frail et al. il99h .1994) , which showed that the H i 
column density varies significantly over scales between 5 and 
110 AU, with 10%-15% of the cold neutral gas distiibuted in 
AU-sized structures (Frail et al. |1994| l. However, some more 
recent radio observations on the same pulsars (Weisberg & 
Stanimirovic 2007 ) have shown that the variations are far smaller 
than those originally found by Frail et al. ( |1991|[T994] l. 

These studies were followed by a series of works looking 
at the variations of Ca ii and/or Na i column densities along the 
lines of sights to close binaries or high proper motion stars (see 
Crawford 120031 and Lauroesch '2007 for a review). Similar in- 
vestigations were carried out for molecular gas (CH, CH^, and 
CN; Pan et al.^004i RoUinde et al. 2003 ), and diffuse interstellar 
bands (Cordiner, Sarre & Fossev 2005) . An alternative method 
is the study of interstellar absorptions along the lines of sight to 
stellar clusters, like M92 (Andrews, Meyer & Lauroesch 1200 11 1 
and Lo-Cen (Van Loon et al. 2 009) . or the Magellanic Clouds 
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(e.g. Andre et al. 120041) . For a general review on the small ISM 
structures in our Galaxy the reader is referred to Haverkorn & 
Goss (I2007]) . 

In this article we present an independent technique to ana- 
lyze extra-galactic ISM structure on spatial scales of about 100 
AU. The proposed method is based on the extremely high expan- 
sion velocity displayed by a supernova (SN) photosphere (^10"^ 
km s"' or 5.7 AU day"'). A Type la SN reaches a photospheric 
radius of ~10'^ cm (~100 AU) in two weeks from the explosion, 
and expands at a rate of 6 to 3 AU day"' during the first two 
months of its evolution. If the typical size of the fluctuations in 
an intervening cloud is much larger than lO'^ cm, then the as- 
sociated absorption features will not evolve with time. On the 
contrary, if the ISM is patchy on comparable scales, the column 
density fluctuations will translate into measurable variations of 
the corresponding absorption features. 

After introducing a simple model for the calculation of time- 
dependent line equivalent widths for a given cloud geometry 
(Sect. j2]i, we present the results of Monte-Carlo simulations 
(Sect. p|, and discuss the applicability of the method and the 
effects on the observations of Type la's (Sect. |4]l. Appendix A 
gives the details on the derivation of the composite equivalent 
width. 



2. A simple model 

Let us fix a reference polar coordinates system {r,6) whose ori- 
gin is located at the explosion center (see Fig.[T]i. As seen from 
a far observer (the distance between the SN and the observer is 
assumed to be much larger than that between the SN and the 



2 



F. Patat et al.: Studying the small scale ISM structure with supemovae 




Fig. 1. Relevant quantities used in the text. The underlying col- 
umn density map was generated using a power-law spectrum 
with y=-2.8 (see Sect.[TT]). 

intervening material), the SN will appear as an expanding pho- 
todisk which extends to the ejecta boundary radius rej(f)- We then 
consider a cloud placed in front of the SN, at a distance large 
enough that the explosion has no effect on its physical condi- 
tions (i.e. >10 pc, see Simon et al. 120091 See also Sect. |4] here), 
and we indicate with N(r,ff) the cloud column density in the 
species under consideration. Finally, we introduce <t)(r, t) as the 
time-dependent surface brightness profile of the photodisk at the 
wavelength of interesj]] This function is normalized as follows 

(D(r, f) rdOdr^l , 

o Jo 

SO that the total continuum flux emitted by the photodisk along 
the line of sight is equal to unity. If ^(A^, b) is the curve of growth 
for the given transition and Doppler parameter b, the equiva- 
lent width {EW) produced by an infinitesimal cloud column with 
cross section dA-rdO dr is dEW - g[N, b] <i>(r, 0) dA. 

If we neglect the small contribution by photons scattered by 
the cloud into the line of sight, the total equivalent width is then 
computed integrating the contribution of each single infinitesi- 
mal element over the photodisk (see Appendix A for the details): 

EW{t) = g[N{r, 0), b(r, 6)] <D(r, f) r dO dr . 

Jo Jo 

So far we have considered the possibility that the Doppler 
parameter b can change across the cloud. However, in the lack 
of evidence for significant variations over the scales of interest 
(see for instance Welty & Fitzpatrick l20011 l. in the following we 
will assume that b is constant across the relevant portion of the 
cloud. We will briefly discuss the efifects of a space-dependent 
Doppler parameter in Sect. [3] 

With the aid of the outlined procedure, one can follow the 
time evolution of EW for any input cloud column density map, 
provided the photodisk's expansion law is known. 

In the assumption of homologous expansion (see for instance 
Jeffery & Branch 1990), the radius of the photosphere rph(f) is 

' We assume 9 has no azimuthal dependence, i.e. that the SN is 
spherically symmetric. For a Type la SN this is a reasonable assumption 
(Wang & Wheeler l2008t . 
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Fig. 2. Surface brightness profiles derived from SYNOW best 
fits of SN 1994D spectra on four different epochs (-10, -4, +7 
and -1-15 days from B maximum). 



obtained multiplying the photospheric velocity Vph(f) by the time 
f elapsed from the explosion. For Vph(0 we adopted a best fit to 
values computed via SYNOW modeling of the spectroscopically 
normal SN 1994D, a standard Type la event (Branch et al. 20051 1. 
During the photospheric phase (f <100 days), Vph(f) is well ap- 
proximated by an exponential law, and the best fit relation takes 
the following form: 

rph(f) = (3.5 + 5.3 e-'/^") f, 

where rph is expressed in AU and t in days from the explosion. 
This law shows that the change in the photodisk dimensions is 
significant during the first two months: from day -10 to day -1-50 
(counted from maximum light) increases by more than a fac- 
tor of 5 (60 AU to 320 AU). 

As for the photodisk surface brightness profile, this was com- 
puted using a modified version of the spectrum synthesis code 
SYNO\\0 (Branch et al. l2005T l. The profile, obtained from best 
fits of SN 1994D spectra for the Nai D rest-frame wavelength, 
rapidly drops at r = rph during the early phases (see Fig. l2]i. 
As time goes by, a significant fraction of the flux (up to ~lo% 
on day -1-15) is emitted above the photosphere, and it is due to 
scattering by the broad Nai D doublet intrinsic to the SN (see 
for instance Jeffery & Branch 1990). At all epochs, <i>-Q for 
r/rph >2.6, which we used as the effective external boundary 
of the ejecta. To include the time dependence we tabulated <l)(r) 
for a number of epochs (-10, -4, +7, +15, +28 and +50 days 
from B maximum) and subsequently used a linear interpolation 
to derive the profile at any given epoch. The time from B maxi- 
mum light was converted into t using the rise time of SN 1994D 
(18 days; Vacca & Leibundgut I 19961 ). In view of the lack of very 
early spectra, we conservatively assumed that <i>-0 for r > rph at 
t=0. 



^ Note that SYNOW, like some other SN spectrum synthesis codes 
(see e.g. Mazzali & Lucv ll993] l. has zero limb darkening. 
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Table 1. Estimated Equivalent Width RMS errors (Tew com- 
puted using Equation[l](5A^/?=100 and A^=10'2 cm^^-. 
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2A. Cloud generation 

We explored two possible cases: i) an isolated, homogeneous, 
spherical cloud with radius rc and offset xc with respect to the 
line of sight; ii) a patchy sheet with an input spectrum for the 
column density fluctuations. 

In the first case, the column density profile is N{r') - 
No \/l ~ ('"'/'"c)^' where r' is the projected distance from the 
cloud center ( r' < rc), and A^o is the column density correspond- 
ing to a ray going through the center of the cloud. The average 
column density is {N}-2/3 Aq. As the case of an isolated, small 
cloud is probably quite unrealistic, it may be rather regarded as 
a simplified model of an over-density on an otherwise homoge- 
neous sheet. 

For the patchy sheet we adopted a procedure similar to 
the one described by Deshpande (2000). Molecular clouds 
(Elmegreen & Falgarone [1996), the diffuse ionized compo- 
nent (Cordes et al. 11991b . and Hi (S tanimirovic et al. 119991 
Deshpande, Dwarakanth & Goss 120001 ) have a fractal structure, 
characterized by a power-law behavior. Therefore, if ^ = Iq/1 is 
the wave-number corresponding to a given spatial scale / (where 
/() is the maximum spatial scale under consideration), the power 
spectrum of the fluctuations can be written as P(k) oc k'>'. For the 
large scale structure of Hi in the SMC y - -3 (Stanimirovic 
et al. 119991 ) and 7=-2.75 for the cold atomic gas in the Galaxy 
(Deshpande et al. 12000). Since the variations observed on scales 
of ~100 AU and below can be explained in terms of a single 
power spectrum description (Deshpande 2000), we computed 
the column density maps using a power-law with y^-l.^. After 
generating a power spectrum map P(k^, ky) in the Fourier plane 
(with k = (kl + k^y^^), we derive the real and imaginary parts 
of the Fourier transform of N{x,y) using '\/P(k) as the modu- 
lus of the complex numbers, while phases are generated as ran- 
dom numbers uniformly distributed between and 27r. The col- 
umn density map is then obtained anti-transforming into ordi- 
nary space. Given the radius of the photosphere during the time 
interval typically covered by observations (<400 AU), we adopt 
;o=1024 AU. 

Once the two-dimensional column density map is generated, 
it is re-normalized to have an average column density (N) and 
peak-to-peak fluctuations AN (with AN < 2(A)). An example 
cloud is presented in Fig. [T] If AA(0 is the column density dif- 
ference between two points separated by a distance /, the numer- 
ical simulations show that for y=-2.8 the RMS value of AN(lo) 



is ~10% of the maximum variation on the same spatial scale. 
Furthermore, the differences at higher wave-numbers decrease 
proportionally to k^-^^^^^^, as predicted by theory (Deshpande 
2000). This implies that the fluctuations expected on scales of 
100 AU are xAQ% of those observed on scales of 10 00 AU , 
which for Na i can reach a few lO'^ cm"^ (Andrews et al. l2001l ). 
Therefore, a significant amount of structure is expected on spa- 
tial scales comparable to the typical photospheric radius. 

2.2. Measurement error estimates 

The ability of detecting small EW variations hinges on the pre- 
cision to which EWs can be measured. In turn, this relates to 
the instrumental setup and the signal-to-noise ratio SNR reached 
on the adjacent continuum per resolution element. To estimate 
the expected RMS uncertainty ctew we have used the following 
formula (Chalabaev & Maiflard [19831 see their Equation A. 10), 
derived for photon noise-dominated conditions and with no as- 
sumption on the line profile: 

where AA is the line integration range, 6A is the dispersion (A 
px"'), = AA/SA is the number of pixels corresponding to 
the integration range, F is the average flux within the absorp- 
tion, Fc is the average continuum level, and cte^. is the associ- 
ated uncertainty. In addition to the parameters directly related 
to the instrumental setup, the application of this analytical ex- 
pression requires the knowledge of other quantities which de- 
pend on the line profile. To compute them for a given column 
density, we generate synthetic Voigt profiles with input velocity 
parameter and a sampling equal to 6A. This is then convolved 
with a Gaussian profile with a given full width half-maximum 
(FWHM) to mimic the instrumental broadening. Using this pro- 
file we calculate the line integration range AA, which we define 
as the wavelength interval in which the equivalent width is 99% 
of total. Within this interval we finally derive the average line 
flux F. As for the continuum F^, we assume this is computed 
on two symmetric, equally extended regions on either side of 
the absorption. If AA^ is the total extension of this region, then 
the number of pixels used to determine Fc is Nc = AAc/6A, so 
that cr^ /Fc - l/{Nc SNR^). In the computations we have used 

AAc=2 A which, given the featureless nature of the SN pseudo- 
continuum, is a realistic valu^ 

The results obtained for ^=10'^ cm 2, and SNR=100 for dif- 
ferent values of FWHM, b (l, 3 and 5 km s"' corresponding 
to EW of 58, 113 and 138 mA respectively ), and 6 A are pre- 
sented in Table [T] As the RMS errors are inversely proportional 
to SNR, these values can be readily scaled to different signal- 
to-noise ratios. These results have been checked against Monte- 
Carlo simulations and were found to be consistent to within a 
few 0. 1 mA. Incidentally, this questions the need for a revision 
of the Chalabaev & Maillard formula discussed by Vollmann & 
Eversberg (2006). 

In the following we will consider an equivalent width vari- 
ation AEW detectable if [AisWI > 5 's/Ictew- For a typical case 
where FWHM^l km s^', km s^', 5/1=0.01 A pix^', and 

^ For the purposes of measuring the EW of an inter-stellar line, the 
continuum definition in a SN spectrum is much easier than in a stellar 
spectrum, where the presence of other intrinsic features may contami- 
nate the adjacent regions. 
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factor, it does not produce a saturated profile in the emerging 
absorption line. 

A distinctive feature of small (rc <64 AU), offset clouds is 
an EW growth followed by a decrease during the pre-maximum 
light phase (Fig. [3] lower panel). This is due to the growth of 
the covering factor as the photodisk starts intersecting the off- 
centered knot. Once a maximum value is reached, the subse- 
quent increase in the photodisk size causes the covering fac- 
tor to drop. Although this mechanism can produce an absorp- 
tion feature which grows in strength and then disappears on 
timescales of a month, this is expected to happen only during 
the pre-maximum epochs, when the relative increase in the pho- 
todisk surface is very fast. Larger clouds placed at larger offsets 
(xc >256 AU) produce features which start appearing around 
maximum light, but keep steadily increasing in strength up to 
several months after maximum. 

Finally, for central column densities (or density contrasts) 
smaller than ^2x10'" cm"-, the EW variations are always below 
the detection limit. 



Fig. 3. Examples of simulated Na i D2 EW variation as a function 
of time for a spherical, homogeneous cloud with offset xc=0 
(upper panel) and xc-6A AU (lower panel). The various curves 
refer to different cloud radii rc (8, 16, 32, 64, 128, 256, 512 
and 1024 AU from bottom to top). In all cases A^o(Nai)=5xlO" 
cm"^ and b- \ km s"'. 

SNR-IQQ, this turns into a 5-cr detection limit A£'W,,„,=4.4 mA 
{AEWii,n^53 mA for ^7=5 km"'). 

3. Results of simulations 

Although the model can be used for any inter-stellar absorption 
line, in the following we present the results obtained for Nai D2, 
because it is a strong transition, it falls in a region almost free 
of telluric absorption features, and in a spectral interval where 
most optical, high-resolution spectrographs have their maximum 
sensitivity. 

3.1. Isolated spherical cloud 

Example EW evolutions for two different cloud offsets xc (0 and 
64 AU) and a number of cloud radii rc are presented in Fig. [3] 
up to 2 months after maximum light. In general, the maximum 
variability is expected when the cloud is close to the center of the 
photodisk. The maximum variation is achieved for cloud radii 
between 64 and 128 AU. Also, small clouds are better detected 
during the early phases (when their size is comparable to that of 
the photosphere), while the detection of large clouds requires a 
larger time span. If the cloud is too large {rc >512 AU), then 
the EW variation is not sufficiently ample to be detected. With 
a minimum set of two observations taken 10 days apart in the 
pre-maximum phases, for typical high-resolution setup, a 5-cr 
detection limit of 4.4 mA, b=\ km s"' and A^o=5xlO" cm"^, 
the simulations show that one is able to detect clouds with rc 
between 16 and 128 AU up to a maximum offset of 64 AU. For 
larger offsets the cloud starts to intersect the photodisk when its 
size is too large and the corresponding covering factor is too 
small. Besides implying probably unrealistic density contrasts, 
increasing the column density does not enhance the detectability 
of a small size cloud, since this rapidly becomes totally opaque. 
We note that, while this causes the saturation of the covering 



3.2. Patchy clouds 

Because of the stochastic geometry of the power-law clouds, 
their effects were evaluated using a statistical approach. For a 
given cloud realization we computed EW{t) and derived the 
absolute peak-to-peak variation A/iW over the whole time in- 
terval -10 to +50 days. To mimic a more realistic situation, 
we did this also on two sub-sets of data-points including only 
post-maximum observations (0, +10, +20, +30, +40, +50 and 
0, +10, +20, +30 days). As statistical estimators we computed 
the median absolute peak-to-peak variation over the full range 
{hEW), the semi-inter-quartile range and the 99-th percentile 
(AiiWgg). For the three time ranges we finally estimated the 5- 
cr detection probabilities P, (-10< t <+50), Pq (f >0) and P30 
(0< t <+30). The results are presented in Table |2] for different 
values of (A^), AA^ and b. For each parameter set 5000 cloud real- 
izations were computed (blank values indicate detection proba- 
bilities <2xl0""^). {EW) indicates the EW coiTesponding to (A^) 
and the input Doppler parameter. 

Visual inspection of a number of realizations showed that the 
EW{t) curves are smooth, with typical timescales of the order 
of 10 days. The variation rate is systematically larger at early 
epochs and decreases as time goes by. The simulations show 
that peak-to-peak column density fluctuations smaller than 10" 
cm"^ on a scale of ~1000 AU do not produce any measurable 
effects. Even with high (~100) signal-to-noise ratio observa- 
tions, starting 10 days before maximum and covering the first 
two months of the SN evolution, the detection probability is be- 
low 2%. This grows significantly when the fluctuations exceed 
~10'^ cm"^. Incidentally, this implies that SNe suffering higher 
extinction are expected to display more pronounced variations 
(see also Chugai l2008l l. Because of the saturation effect, the EW 
variations are more marked for larger values of the Doppler pa- 
rameter. Finally, for column densities exceeding ~10'^ cm the 
detection probability decreases due to line saturation. 

To study the effect of a spatially variable Doppler parame- 
ter, we have run a set of simulations in which b is allowed to 
fluctuate around the average value across the cloud. We have 
tentatively modeled the Doppler parameter map using the same 
algorithm and spatial scales spectrum adopted for the column 
density generation (see Sect. |2.1[ i. We remark that this is not 
meant to reproduce real physical conditions, but only to esti- 
mate the consequences of velocity dispersion fluctuations (for 
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Table 2. Results of Monte-Carlo simulations for power-law 
clouds with y=-2.8. Equivalent width variations are peak-to- 
peak (see Sect. |3. 2 1. 
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instance, b and (A^) were left completely independent). The MC 
runs (with 1< b <5 km s"') show that for (A^) <10" cm"^ there 
is practically no difference with respect to the case with con- 
stant Doppler parameter. The differences start to be significant 
at (A^) >10'^ cm"-, for which the typical variations and detec- 
tion probabilities become much larger. 

This was to be expected, since in the quasi-linear regime at- 
tained at low column densities the effect tends to average out. 
On the contrary, as one enters the non-linear part of the curve of 
growth, the largest variations are produced by the regions of the 
cloud where the Doppler parameter is higher (i.e. less subject to 
saturation), thus skewing the distribution towards more marked 
EW fluctuations. The exact behavior depends on the way the ve- 
locity dispersion varies across the cloud and how this (if any) 
relates to the column density fluctuations. However, the conclu- 
sion that a variable Doppler parameter enhances the detection 
probability is of general validity. 



4. Discussion and Conclusions 

The simulations presented in this paper indicate that marked 
time effects on the measured EWs are expected for small 
(rc <100 AU), isolated clouds with M)(Nai) -lO^-lO'^ cm^^ 
and for small offsets {x^ <100 AU). However, the existence of 
such structures is seriously questioned in terms of pressure equi- 
librium arguments and the yet unknown processes that would 
produce them (see Heiles [19971 and references therein). Frail et 
al. (1994) detected maximum AA^(H i) variations that range from 
~10''' to ~5xlO^° cm"^ on scales between 5 and 100 AU. For a 
Galactic Nai/Hi ratio (Ferlet, Vidal-Madjar & Grv |1985l l fliis 
turns into AA^(Nai) between 5xlO'° and 3x10'^ cm"^. These 
large changes have been interpreted as arising within ubiqui- 
tously distributed small structures. However, this picture has 
been questioned by Deshpande (2000), who has convincingly 
shown that the observations are consistent with a single power- 
law description of the ISM, down to AU scales. In these cir- 
cumstances, peak-to-peak variations AA^(Nai) ~10'^ cm"^ are 
expected on scales of ~100 AU. Our calculations (see Table |2]l 
show that for a Type la SN observed under the most favorable 
conditions {SNR>10Q on all epochs, spanning from -10 to +50 
days) these would appear in less than 10% of the cases for b ~l 



km s"'. This fraction increases to ~80% for b ~5 km s" , but 
in all cases it is AEW <40 mA. These small variations imply 
negligible changes in E(B - V), but can have some effect on 
observing programmes studying the evolution of Nai features 
possibly arising in the circumstellar environment of Type la pro- 
genitors (Patat et al. 120071 Simon et al. l2009l l. However, we note 
that these values are more than a factor of 10 smaller than the 
Nai D variations detected in the Type la SN 2006X (Patat et al. 
[2007 ), which were attributed to the ionization effects induced by 
the SN on its circum-stellar environment. In contrast, no statisti- 
cally significant variations were detected for the CN, CH, CH^, 
Ca I lines and DIBs associated to an interstellar cloud in the host 
galaxy (Patat et al. 2007. Cox & Patat 2008). In the only other 
well studied case published so far (SN 20071e), the EW of four 
Na I D components remained constant to within a few mA dur- 
ing six epochs spanning about 3 months (Simon et al. 20091). 
In this time range rph changed approximately from 100 to 400 
AU, and the lack of evolution is in line with the predictions of 
our model for a power spectrum ISM and definitely excludes 
the presence of small, isolated clouds with sizes comparable to 
rph. Although the available data are still scanty, the multi-epoch, 
high-resolution campaigns which are being conducted for the 
study of Type la progenitors will provide a more statistically sig- 
nificant sample. 

All the discussion so far is based on the assumption that the 
physical conditions of the ISM are not modified by the SN explo- 
sion, so that all variations in the absorption lines are due to pure 
geometric effects. Indeed, a Type la SN can produce changes 
in the ionization balance of low-ionization species (like Na i or 
Ki) up to quite large distances, of the order of 10 pc (Patat et 
al. 2007, Chugai ,2008, Simon et al. 2009). Giv en the UV flux 
predicted for these distances (Simon et al. I2009I I, the ionization 
timescale of Nai is expected to be about 120 days. Since the 
electron density in the ISM is low, the recombination time is 
extremely long. Therefore, under the assumption of a constant 
ionizing field, one would expect the amount of neutral Na to de- 
crease with timescales of several months, hence mimiking small 
scale structure effects. However, the UV flux of a Type la SN 
decreases significantly after maximum light (a factor ~20 in the 
first 40 days; Brown et al. 2009) implying that the maximum 
distance is probably less than 10 pc. Given the range of possi- 
ble distances to an inter-stellar cloud within the host galaxy, this 
suggests that significant Nai column density variations in the 
ISM induced by the SN radiation field are improbable. 

Another important fact is that, because of the much higher 
ionization potential of Can, and the strong UV line blocking 
present in Type la spectra, the EW of the ubiquitous H&K lines 
becomes insensitive to the SN radiation field already at distances 
of a few 0.1 pc (Simon et al. 2009). In contrast, in the case of a 
geometrical origin, all species are expected to show synchronous 
variations, although possibly with different amplitudes and time 
scales (Lauroesch & Meyer 2003 ). Therefore, a comparison be- 
tween the behaviors of Nai and Can should allow one to dis- 
entangle between geometrical and ionization effects, similar to 
what has been proposed by Patat et al. ( 120071 ) . 

The method we presented enables the study of small scale 
structure in the extragalactic ISM, which is out of reach for any 
of the techniques deployed so far. In this respect we note that 
the same method can be in principle applied to the Galactic ISM 
to study its sub-AU structure, for which no direct measurements 
are available yet. Although probably requiring very high signal- 
to-noise ratios, this technique might put important constraints on 
the very small scale structure. In this article we have discussed 
the case of the strong, easily detectable Nai D lines. However, 
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Other weaker lines (e.g., Ki, Cai, Can) can be used, especially 
when the Na i D lines are saturated. In general, the simultaneous 
study of different atomic/diatomic lines along the lines of sight to 
SNe will contribute to get a more detailed picture of the physical 
conditions of the ISM in the small scales regime. 

Appendix A: Calculation of composite equivalent 
width 

Let us consider an extended photodisk, a cloud placed in front of 
it and an absorption line with profile function (p{v), normalized 
so that 



f 

Jo 



(v)dv^ 1. 



If is the space-dependent column density of the cloud 
(atoms cm"^) for the atomic species under consideration, the 
monochromatic line optical depth is t(v) - K N <piv), where K is 
the frequency-integrated line opacity (cm^ atom"'). Let us now 
consider an infinitesimal element of the photodisk dA - r dO dr 
in the polar reference system introduced in Sect. |2] (see Fig. [T]i, 
and define <b{r) as the radial profile of continuum surface bright- 
ness for the photodisk. For the sake of simplicity we assume that 
the continuum Iq is constant across the line profile and is nor- 
malized to unity: 



/o = 27r r ' (5(r) r d6 dr ^ \ . 
Jo 

With these settings, the monochromatic intensity contributed 
by the infinitesimal photodisk element dA to the total observed 
intensity can be written as 



dl{v, r, 0) = e 



-K N(r,8) 0(v) 



(D(r) rdddr . 



For a distant observer, to whom the photodisk will appear 
as an unresolved source, the total line intensity profile I{v) is 
obtained integrating all the infinitesimal contributions over the 
photodisk: 



/(v)= e 
Jo Jo 



-K N(r,9) ipiv) 



<E)(r) rd0dr. 



EW 



Given the definition of equivalent width 



Jo 



1 - 



dv , 



the composite equivalent width can be expressed as follows: 

'S>ir)rdedr\ dv . 



Jo \ Jo Jo 



which implies that the composite equivalent width is the 
weighted sum of the equivalent widths produced within each in- 
finitesimal cloud element. If ^(A^, b) is the curve of growth for 
the given transition (where b is the Doppler parameter that char- 
acterizes the line profile), then Equation |A. 1 [ can be reformulated 
as: 



Because of the normalization of <l'(r), the previous relation 
can be rewritten as 

p/Vj p2:T p+oo 

EW ^ \ [l-e-'^'^^'''^^'^^'^) mr,0)rdedrdv{A.l) 

Jo Jo Jo 

Now, the inner integral is the equivalent width one would 
observe if the physical system were composed only by the in- 
finitesimal cloud element dA and the corresponding portion of 
the photodisk: 



EW{r, 6) 



■ e 



-K N(r,0) <p{ 



dv . 



EW 



- f" 

Jo Jo 



g{N(r, ff), b{r, 0)] 0(r, ff) r dO dr 



(A.2) 
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